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We develop a gauge-invariant formalism to describe metric perturbations in five-dimensional 
brane-world theories. In particular, this formalism applies to models originating from heterotic M- 
theory. We introduce a generalized longitudinal gauge for scalar perturbations. As an application, 
we discuss some aspects of the evolution of fluctuations on the brane. Moreover, we show how the 
five-dimensional formalism can be matched to the known four-dimensional one in the limit where 
an effective four- dimensional description is appropriate. 



I. INTRODUCTION 
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In recent years, the way in which string-theory is believed to be connected to observable physics has changed 
ON ' dramatically. The new viewpoint is mainly due to two ideas, namely the brane-world idea and the idea that 

couplings and scales of additional dimensions are much more flexible than previously assumed |^,^,^U|]. Not only 
have these ideas led to new directions in M-theory phenomenology and, more generally, string-theory inspired particle 
phenomenology, but also in early universe cosmology. 
. Much of the recent activity in brane-world cosmology is centered around five-dimensional brane-world theories, see 
' for example ||-|37|. A large class of such theories arises from heterotic M-theory Q^,^]. Other five-dimensional 



models have been introduced in ref. [ff0||-|43[| which may provide an alternative solution to the hierarchy problem. 

A central question is whether the possible existence of a brane-world and large additional dimensions in the early 
universe leads to observable consequences today. Specifically, cosmological perturbations as, for example, observed 
in the cosmic microwave background provide us with a window to the early universe that, perhaps, can be used 
to test the brane-world idea. It is with this motivation in mind, that we set out to study metric perturbations 
in brane-world models. It may not be immediately clear that the existence of additional dimensions and branes 
should have important consequences for the formation and evolution of cosmological perturbations. Let us, as a 
[ comparison, consider "traditional" Kaluza-Klein cosmologies, where the higher-dimensional universe is usually split 
i into a product of maximally symmetric subspaces each one with an individual scale factor. Cosmological perturbations 
are normally treated at the linearized level and, hence, in such Kaluza-Klein cosmologies [ff7|-|4^] internal and external 
• i-h , fluctuations basically decouple. The additional dimensions effect perturbations of the three-dimensional universe only 
via the kinematics of the background scale factors (and/or dilaton background fields @). The situation is quite 
different for brane-world theories. The branes are stretched across the three-dimensional universe and are located 
at specific points in the internal space. Moreover, they carry worldvolume fields that can only propagate on the 
brane and that are likely to be excited in the early universe, both coherently and thermally. As a consequence, the 
cosmological background in an early brane-world universe is highly inhomogeneous in the additional dimensions since 
the branes constitute localized sources of stress-energy. Even if perturbations around such a background are treated 
at the linearized level, the perturbations of the three-dimensional universe are effected by the non-linear distortion of 
the cosmological background in the internal dimensions. This constitutes a crucial difference between conventional 
Kaluza-Klein cosmology and brane-world cosmology which is directly related to the presence of branes. It is this 
difference that may lead to new predictions for cosmological perturbations in brane-world models and that motivates 
the present investigation. 
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In this paper, we will develop a formalism for metric fluctuations in brane-world theories that takes the charac- 
teristic property of brane-world cosmologies, the above mentioned inhomogeneity in the additional dimensions, into 
account. We understand such a formal development as a first step towards analyzing predictions for cosmological 
perturbations in brane-world theories. In the next section, we start out to generalize the well-known formalism of 
four-dimensional gauge-invariant metric perturbations E3] , p5[ to brane-world theories with an arbitrary number of 
additional dimensions. Subsequently, in section III, we focus on five-dimensional brane-world models on the orbifold 
S 1 jZi related to those originating from heterotic M-theory jl 38 3jJ. Specifically, we consider the five-dimensional 
Einstein equation coupled to bulk as well as brane stress-energy. For such a theory we derive the equations of motion 
for scalar metric perturbations in a generalized longitudinal gauge. Consistency of these equations is used to determine 
the most general structure of the stress-energy on the brane. In section IV, these results are applied to find the Israel 
matching conditions Q] for the scalar metric perturbations restricted to the branes. Finally, in section V, we show 
how our formalism for brane-world metric perturbations is related to the conventional one in four-dimensions. This 
is done by a matching procedure applied in a limit where the five-dimensional brane-world theory has an effective 
four-dimensional description. 



II. GAUGE-INVARIANT VARIABLES 



In this section, we will develop a gauge-invariant formalism for metric perturbations in brane-world models. Using 
such a gauge-invariant approach is particularly useful in order to identify the correct physical degrees of freedom. 
Once this has been done, a specific gauge can be chosen in order to simplify the subsequent equations. Specifically, 
we will use a generalized longitudinal gauge later on. However, as a warm-up for the higher-dimensional case, we 
would first like to review the well-known four-dimensional gauge- invariant formalism following ref . [fi"5| . 



A. The four-dimensional formalism 



Starting point is a background metric with a maximally symmetric three-dimensional spatial space. This metric is 
of the form 

ds\ = a\ {dt 2 - <>,,</r',/r<} , (1) 

where indices i, j, ■ ■ ■ — 1,2,3 run over the three spatial indices. Indices /i, v, ■ ■ ■ = 0, 1, 2, 3 are used to index four- 
dimensional space-time. Furthermore, a 4 = a 4 (t) is the four-dimensional scale factor and f2y is the metric of the 
three-dimensional maximally symmetric space explicitly given by 

-2 

(2) 

Here k = 0, 1, —1 corresponds to a flat, closed or open universe, respectively. The idea is now to classify perturbations 
of the metric (|l|) according to their transformation properties with respect to the maximally symmetric space. This 
leads to the perturbed metric 

dsj = a\ {(1 + 2<j) 4 )dt 2 - [(1 - 2^ 4 )fiij + 2£ 4 |y + 2F 4m + h 4lj ] dx l dx ] + W 4l dtdx 1 } . (3) 

Here and in the following four-dimensional quantities are indexed by "4" to distinguish them from their higher- 
dimensional counterparts to be introduced later. The vertical bar refers to a covariant derivative with respect to 
the metric fly. The vector F 4 i has a vanishing divergence, that is F 4 j} 1 = and the tensor haj is traceless and 
divergence- less, that is ha = and h 4i y = 0. In addition, we can decompose the off-diagonal perturbation Wu 
further into the gradient of a scalar B 4 and a divergence-less vector Sa . Explicitly, this reads 

W 4l = B 4[i + S 4l . (4) 

Consequently, we have four scalar metric perturbations (</>4, ip 4 , E 4 , B 4 ), two vector perturbations (F 4 i,S 4 i) and a 
tensor perturbation (h 4i j ) . All these perturbations are functions of time as well as of the spatial coordinates x l , of 
course. Next we consider an infinitesimal coordinate transformation 

x" -> = x" + ^ , (5) 



I + -kx l x m S h 
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where the vector £ M depends on all four coordinates, in general. The corresponding infinitesimal change of the metric 
is given by 



<?4a«/ — * 94,/mu = g±fiu — 2V( M £„) . (6) 

To understand how this coordinate transformation acts on the metric perturbations we split £ M as £ M = (£°,£; 4 ) into 
a time and a spatial part. The spatial component can be decomposed further into a gradient and a divergence- less 
part as 

C=^ + r, i . (7) 

As a result, the transformation parameter £ p contains two scalar components (£ , £) and one vector component (C)- 
Given this setup, one can compute the transformation properties of the metric perturbations by applying eq. (|6|) to 
the perturbed metric (^|) and taking into account that £p = a 4 (£o, — For the scalar perturbations one finds 

04 = - H 4 e - i° , (8) 

V^ 4 = ^4 + , (9) 
B 4 = B 4 + C° - i , (10) 

£ 4 = ^4 ~ £ ■ (11) 

Here, -H4 is the Hubble parameter defined by H 4 = 04/04. The vector perturbations transform as 

F 4l = F 4l - m , (12) 
S 4l = Su - f]i , (13) 

while the tensor perturbation haj is invariant. In these equations, spatial indices are lowered and raised with the 
metric flij, that is, for example £u = flij^ and rji = Qijrf . With these results, it is straightforward to introduce the 
following gauge-invariant variables. 

Scalar variables 

$4 = <p 4 + H 4 (B 4 - £4) + B 4 - E 4 (14) 
* 4 = ^4 - H A (B A - E 4 ). (15) 



Vector variables 

Ti = S 4i - F 4l . (16) 

Tensor variables 

huj ■ (17) 

As the physical degrees of freedom, one has therefore identified two scalar perturbations, one vector perturbation 
and one tensor perturbation. Of par ticular importance are the two scalar perturbations $4 and ^4 on which we are 
going to focus. The expressions ( jbil) and ( |l5| ) for these perturbations suggest the gauge choice B 4 = E 4 — in the 
scalar sector which is referred to as longitudinal gauge. Clearly, from the above transformation properties of the scalar 
perturbations such a choice can be made. Then, the gauge-invariant scalar variables coincide with the "original" 
variables, that is, $4 = <p 4 and ^4 = -04- This gauge choice considerably simplifies subsequent calculations and its 
generalization will be quite helpful to deal with the higher-dimensional case. The perturbed metric then takes the 
form 

ds\ = a\ {(1 + 2(j) 4 )dt 2 - (1 - 2ijj 4 )n ij dx l dx j } . (18) 

Finally, we need to specify the stress-energy. For the background, by the maximal symmetry of the three-dimensional 
spatial space, it is dictated to be of the form 

T 4 M „ = diag(/9 4 , -P4, -P4,, ~Pi) ■ (19) 
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The scalar perturbations to this stress-energy can be written as 

m »„=( J" „ ^ + /i a ^) (20) 
V (Pi +Pi)a i v^ 1 -dp4<5 j + 0-4' |j / 

with the potential V4 for the velocity field and the quantity a specifying the anisotropic stress. The equations of 
motion for the background and the scalar perturbations subject to the above stress-energy are given in ref. p5[ and 
will not be repeated here. These equations form the basis for the study of cosmological perturbations and we now 
turn to develop their higher-dimensional generalization. 

B. Gauge-invariant variables in brane- world theories 

We would now like to proceed in close analogy with the four-dimensional case reviewed above and develop a 
gauge- invariant formalism of metric perturbations in brane- world theories. First, we consider the general situation 
of d additional dimensions although later we will be more specific and focus on the case d = 1, that is, a five- 
dimensional universe. The coordinates of the additional dimensions are denoted by (y 5 ,- • • ,y 4+d ). For the purpose 
of this subsection, all we need to specify is that the branes are stretched across the usual four-dimensional space-time 
and are located at specific points (or submanifolds) in the additional dimensions. We will be more precise about the 
brane positions later when we consider the five-dimensional case. 

How should metric perturbations be classified in such a brane-world theory? In the previous four-dimensional case 
we have used their tensor properties with respect to the three-dimensional spatial subspace for this classification. At 
first glance, one might now want to use their tensor properties with respect to the 3 + rf-dimensional spatial space. 
The cosmological principle, of course, only asserts the maximal symmetry of the usual three-dimensional space but 
the maximal symmetry of the e?-dimensional internal space may be taken as an additional, simplifying assumption. 
It is at this point, that the brane-world nature of the theory comes into the game. Since the branes are localized in 
the additional dimensions the assumption of maximal symmetry cannot, in general, be extended to those dimensions. 
In fact, as will become more explicit below, the branes lead to stress-energy in the Einstein equation localized in the 
additional dimensions and, hence, the symmetry of the background metric will typically not be enhanced with respect 
to the four-dimensional case. Consequently, we split the coordinates into two groups, namely the inhomogencous 
coordinates (y a ) — (i, y 5 , ■ ■ ■ , y 4+d ) on which the background metric generally depends in a non-trivial way and the 
usual three spatial coordinates (x l ) corresponding to the maximally symmetric space. In the following we use indices 
a, b, ■ ■ ■ = 0, 5, • • • , 4 + d for time and the additional dimensions, indices = 1,2,3 for the three-dimensional 

space and indices a, (3, - ■ ■ = 0, 1, 2, 3, 5, • • • ,4 + d for the full 4 + d-dimensional space-time. Then the most general 
higher-dimensional metric consistent with the maximally symmetric three-dimensional spatial manifold is given by 

ds 2 = a 2 {-f ab dy a dy b ~ <>,,,/.,•',/.,•• } , (21) 

where the scale factor a and the metric 7 a {, are functions of the coordinates y a only. Here fij,- is the metric of the 
maximally symmetric space of constant curvature given in (|^). Given this structure of the background metric, we are 
forced to classify metric perturbations by their three-dimensional tensor properties as in the four-dimensional case. 
We stress again that this is a direct consequence of the brane-world nature of the theory that we are considering. 
With these remarks in mind, the higher-dimensional generalization of the perturbed metric (|j|) can be written in the 
form 

ds 2 = a 2 { lac (5 c b + 20g) dy a dy b - [(1 - 2i/>) fig + 2E {ij + 2F m + ha] dx'dx 1 - 2W ai dy a dx i } . (22) 

As in the four-dimensional case, Fi and hij have a vanishing divergence and, in addition, hij is traceless. As before, 
the three-vectors W a i can be split as follows, 

W ai = B a[i + S ai , (23) 

where SVij = 0. Observe that the perturbed metric (p2|), defined in this way, is completely general. In fact, this 
can be easily seen by counting degrees of freedom. As an example, we can consider the simplest case of only one 
extra dimension setting y = y 5 and a,b, ■ ■ ■ = 0,5. Then, the most general perturbed metric contains 15 degrees of 
freedom, which are parameterized by the seven scalar perturbations ((/>$, 0°, 0|, ip, E, B a ), six degrees of freedom 
from the vector perturbations S a %) and two degrees of freedom from the tensor perturbation hij. Of course, in 
counting the degrees of freedom originating from vector and tensor perturbations we have taken the constraints on 
these quantities into account. 
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Let us now return to the general case of d additional dimensions and consider the coordinate transformations 

x a _> x a =x a +C (24) 

with infinitesimal parameters £ Q . In accordance with the above discussion, we split these parameters as = (£ a , 
We adopt the useful convention that indices of type a (i) are lowered, raised and contracted using the metric ^ a b 
(f2ij). Furthermore, we take the vertical bar to denote the covariant derivative with respect to j a b or f2y depending 
on the index type. From the transformation law 

g a /3 =9ap- 2V (25) 

of the metric, and taking into account that £ a = a 2 (£ a , — we find for the transformation of the scalar perturbations 

$0ab = -f(o|6) ~ -ff C Cc7ah, (26) 

H = H a £ a , (27) 
5E = (28) 
SB a = £ a - £|„. (29) 

where we have introduced the generalized Hubble parameters 

a \c , s 

ff c = — • 30 

a 

The vector perturbations in the metric (122) change according to 



SFi = -77,, (31) 
8S ai = -»7i| a . (32) 



Finally, the tensor perturbation is invariant under the first order gauge transformation (25). From these results 
we easily find the following gauge-invariant quantities. 

Scalar variables 

= <t>ab + H C {B C - E\ c ) lab + (B (a - E l(a ) ]b) (33) 
* = V - H c (B c - B, c ) . (34) 



Vector variables 

Fai = S ai - Fi\ a . (35) 

Tensor variables 

hij . (36) 

We conclude that the physical degrees of freedom consist of the (d + l)(d + 2)/2 + 1 gauge invariant scalar pertur- 
bations ($a&) ^ / ) ; ^ + 1 gauge invariant vector perturbations T a i and a gauge invariant tensor perturbation hij. The 
above gauge-invariant variables are a direct generalization of the corresponding four-dimensional ones. Specifically, 
restricting to no additional dimensions and setting 700 = 1, eqs. (|33|)-(p6|) exactly reproduce the four-dimensional 
expressions (|i^-(|l7|). However, in the case d > our formalism clearly has a richer structure than the conventional 
four-dimensional one. 



C. A generalized longitudinal gauge for scalar perturbations 

In the subsequent sections we will focus on the evolution of scalar perturbations. Vector and tensor perturbations 
will be discussed elsewhere. In order to simplify this discussion we introduce a generalized longitudinal gauge for the 
scalar perturbations. In analogy with the four-dimensional case, this gauge is specified by 

B a = , E = 0. (37) 
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Setting these quantities to zero can indeed be achieved by an appropriate choice of the scalar transformation parameters 
£ a and £ in the eqs. (|2^)-(|29|). Note that we have exactly the correct number of transformation parameters to do 
this and that, consequently, the gauge ambiguity in the scalar sector is complete eliminated by this choice. Then, the 
scalar part of the metric takes the simple form 

ds 2 = a 2 {"fac(5b + 2(j) c b )dy a dy b - (1 - 2tp)n ij dx i dx j } . (38) 

Furthermore, in this gauge, the scalar perturbations <p a b and ip coincide with their gauge-invariant counterparts, that 
is 

$ab = <t>ab, * = 1> (39) 

as it is the case in four dimensions. 



D. The five-dimensional case 



Let us restrict in this section and for the rest of the paper to the case of a single extra dimension y — y 5 . Then, 
the indices a, b, ■ ■ ■ run over the values 0, 5 only. Furthermore, in order to be more explicit, we choose the conformal 
gauge 



( 7ot ) = 6 2 diag(l,-l) 



(40) 



for the background metric j a b by performing a large gauge transformation. Here b = b(t, y) is a new, independent 
scale factor. Then, the perturbed five-dimensional metric (p2|) reduces to 



ds 2 = a 2 {b 2 [(1 + 2(j))dt 2 - 2Wdtdy - (1 - 2T)dy 2 } 

- \p,ij{l - 2-0) + 2E\ij + 2F m + hij] dx l dx j - 2W oi dtdx i - 2W 5i dydx i } , 



where we have defined 



*8, 



*5, 



W — 20q = -2$ 



(41) 



(42) 



Recall that the scale factors a and b are functions of the coordinates t and y only while the perturbations depend on 
all spacetime coordinates. The scalar gauge-invariant variables defined in eq. ( |33| ) can now be written more explicitly 
as 



1 



*i = $° = + - (Ho- Ho)(B -E)+B -E-(H 5 + H 5 )(B 5 - E') 



* 2 = <&f = r- s 



$3 = $0 = 



w 



1 



2 2b 2 

$ 4 = * = i> - y> r 



{H 5 - H 5 )(B 5 - E') + B' 5 - E" - (H + n )(B - E) 

(B'q + B 5 )-2E'- 2H 5 (B — E) — 2H (B 5 - E' 
H (B - E) - H 5 (B 5 - E' 



(43) 
(44) 
(45) 
(46) 



Here and in the following the dot (prime) denotes the derivative with respect to time (the coordinate y). Furthermore, 
we have introduced a second set of "Hubble" -parameters TL a = b\ a /b. Let us specialize these results to the generalized 
longitudinal gauge defined by Bq — B§ — E = 0. Then the above scalar gauge-invariant variables coincide with </>, "0, 
r and W . The metric simplifies to 



ds 2 



a 2 {b 2 [(1 + 2cj))dt 2 - 2Wdydt - (1 - 2T)dy 2 } - (1 - 2^)n i:j dx l dx j } 



(47) 



This metric will be the starting point for our treatment of scalar perturbations in the following sections. In addition to 
the perturbations <f> and ip that we are familiar with from the four-dimensional case it contains two new perturbations, 
r and W, that are related to the presence of the fifth dimension. 
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III. THE PERTURBED EINSTEIN EQUATION IN THE LONGITUDINAL GAUGE 



As we have previously mentioned, the main application we have in mind for this paper is a compactification of a 
five-dimensional theory on the orbifold S\jZ 2 . We start by compactifying the fifth dimension on a circle restricting the 
corresponding coordinate y to the range y G [— -R, R] with the endpoints identified. The action of the Z 2 orbifolding 
symmetry on the circle is taken to be y — > —y. Consequently, there exist two fix points at y = y\ = and y = 2/2 = R- 
We assume that the three-branes, stretching across 3 + 1-dimensional space-time, are located at these fix points in the 
orbifold direction. This setup is appropriate for a large class of five-dimensional heterotic M-theory models fl ^8|j39| 
that originate from 11 -dimensional Hofava-Witten theory. It also applies to the five-dimensional models introduced 
in ref. 

Next, we should truncate the five-dimensional metric in order to make it consistent with the orbifolding. Since the 
metric has to be intrinsically even under the Z 2 action its various components satisfy the constraints 

9nv{-y) = 9i*u(y), (48) 

g^(-y) = -.9^5 (y), (49) 
.955 (-y) = 555(2/)- (50) 

At the same time, we have to make sure that coordinate transformations do not lead out of the class of metrics defined 
this way. The parameter £ Q for an infinitesimal coordinate transformation has, therefore, to be restricted by 

e(-y) = e(y), (si) 
ZH-y) = -Z 5 (y) , (52) 

which directly follows from (p5|). From these rules we can deduce the Z2 properties of the various quantities in 
metric (^) for scalar perturbations. While the background scale factors a, b as well as the perturbations <j>, ip and T 
are Z 2 even, that is, for example, a{—y) = a(y), the perturbation W is Z 2 odd, that is W(— y) = —W(y). Similarly, for 
the scalar components in the transformation parameter £", we find that £0 and £ are even while £5 is odd. Also note 
that the derivative along the fifth dimension of an odd variable is even and vice versa. For instance, W'(y) = W'(—y). 
Based on these Z 2 truncations we should now discuss the continuity properties of all quantities. Normally, one 
requires the metric to be continuous across the whole of space-time in order to have a sensible notion of length and 
time. We will also adopt this viewpoint, however with an additional subtlety. Since the orbifolding identifies the 
upper and lower half of the circle in the fifth dimension already one of them, say the upper half, constitutes the whole 
of space-time. In fact, instead of working with the orbifold picture where one keeps the full circle as we do here, one 
could also use the boundary picture where only one half of the circle (a line-segment) is considered. This shows that a 
jump of a metric component at an orbifold fix point does not contradict the continuity requirement. Of course, such 
a jump is possible only for an odd component of the metric. Concretely, we therefore require that all components 
of the metric ( fl7j ) are continuous across the full orbicircle except for the odd component W which may jump at the 
fix points (but is continuous otherwise). Corresponding assumptions have to be made for the parameter £" so that 
coordinate transformations do not change these continuity properties of the metric. Clearly the even components £ 
and £ have to be continuous then. Is the odd component £5 allowed to jump at the orbifold points? Eq. (|2^) shows 
that r = —05, which has to be even and continuous, transforms with the derivative £ 5 . Hence, if £5 jumped at the fix 
points it would lead to a delta-function singularity in the metric which is clearly unacceptable. We, therefore, have 
to require that £5 is continuous everywhere on the orbicircleQ In particular, this means that £5 vanishes at the fix 
points, that is £5(2/71) = 0. 

It is clear that the above conclusions depend somewhat on the fact that we are working with an orbifold. For 
example, if we had considered compactification on a circle instead, all components of the metric had to be continuous. 
Correspondingly, some of the conclusions below will be slightly modified for other compactifications, however, in a 
way that is usually rather obvious. 

Given this setup the Einstein equations can be written as 

1 2 

G at 3 = R a /3 - ^9apR = T af3 + ^2 T !£l S (y ~ Vn), (53) 

71=1 



1 From this conclusion we see, that we have glossed over a subtlety when introducing the generalized longitudinal gauge. 
Clearly, for continous £5 the quantity B5 can only be gauged to zero if E\g — B5 is continous, as can be seen from ( Bq ) and 
(p9|). We will, therefore, in addition require the continuity of E\ 5 — B$. 
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where we have set the five-dimensional Newton constant to one, for simplicity. The delta-functions in this equation are 
covariant with respect to the fifth dimension, that is, they include a factor of 955. Furthermore, T a p is the bulk 
stress-energy tensor induced by fields that propagate in the full five-dimensional space time. The brane stress-energy 
tensors , on the other hand, originate from fields that are confined to the branes at the orbifold fix points. 

In order to proceed further, we need to specify these stress-energy tensors. Two requirements should be taken 
into account when doing this. First, one should use the fact that the background has a maximally symmetric three- 
dimensional space. Secondly, the brane stress-energy tensors should be restricted in a way that is consistent with 
their nature of representing fields on the branes. This latter requirement can be most easily implemented by using 
the constraints that follow from the Einstein equation ([53]) itself. Concretely, the delta-functions on the right-hand 
side of this equation have to be matched by corresponding delta-functions that appear on the left-hand side. The 
appearance of these latter delta-functions, however, is controlled by the structure of the equations and the continuity 
assumptions about the metric discussed above. 

Let us see what this implies in detail. We start with the background stress-energy. For the bulk, the most general 
form of this tensor consistent with the three-dimensional maximal symmetry is 





T a fj — -pFj . (54) 
\r -qj 

In particular, we note that the 05 component can be non- vanishing. This possibility is, in fact, already realized for 
the simple case of a bulk scalar field that depends on t and y. As far as the symmetry of the background metric is 
concerned, the background brane stress-energy tensors should have the same structure as (|54|). However, as we will 
see in a moment, there are two more requirements that follow from the equations of motion, namely that the 55 and 
the 05 components vanish. As a result, the background stress-energy on the branes has the form 

p (n) 

T ; = ( -p^S'j J . (r,r» 

Let us now proceed to the perturbed stress-energy tensors. Since we are focusing on scalar perturbations we write 
the most general perturbation of the background bulk tensor (jU) that can be expressed in terms of scalars on the 
maximally symmetric subspace. This leads to 

Sp — (p + p)b~ 2 v\j — Sr \ 

ST a p =\ (p + pW* -Sptfj+a^ -uM (56) 

Sr + 2r{<j> + T) - (p + q)W -b~ 2 u\j -Sq J 

where v and u are two potentials for "velocity" fields and a, satisfying cr^u = 0, determines the anisotropic stress. The 
perturbed brane stress-energy tensors should have the same structure. However, as we will see below, the equations 
of motion impose further constraints implying vanishing 55 and 5z components as well as vanishing anisotropic stress. 
Therefore, the brane stress-energy perturbations are given by 



5T {n)a , 



1 £p(") -(p(«) + p ("))b- 2 v\f -<JrM 

-SpM&j 3 I (57) 
^ Sr^-p^W 



We would like to present the equations of motion based on the metric ( |47| ) and on the above stress-energy tensors 
that follow from the Einstein equation (|53|) . However, in writing the metric (47) two gauge choices where involved 
and it is not, a priori, clear that these choices can be made while, at the same time, keeping the branes at y = const 
hypersurfaces as we have conveniently assumed in writing eq. (|53|). As the first choice, we decided to work with 
the two-dimensional metric "f a b in conformal gauge. Fortunately, this can be achieved while keeping the branes at 
hypersurfaces y = const |fL3| . In addition, we have chosen the generalized longitudinal gauge ( |37j ) for the scalar 
perturbations. A brane, described by y = y n before the gauge transformation that leads to longitudinal gauge, is 
described by y = y — £ 5 (y) = y n after this gauge transformation, where y is the transformed y coordinate. However, 
since £ 5 (j/n) = 0, as discussed above, this equation is solved by y = y n to linear order. We conclude that, in the new 
coordinates that lead to the generalized longitudinal gauge, the brane location is unchanged to the relevant linear 
order. In summary, therefore, using the Einstein equation ( |53"| ) with the branes described by y — y n does not restrict 
the generality of our results given the gauge choices that we have made. 
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The background equations following from (53) have been given in ref. |13] for the case of stress-energy induced 
by scalar fields and in ref. JTHJ for the case of ideal fluids. For completeness and in order to incorporate some of 
the generalizations that we have made (such as the inclusion of a 05 component of the bulk stress-energy) we will 
nevertheless present these equations here. We find 



a 2 b 2 G° = 3 



aVG 5 5 = 3 



ab 
ab 
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(58) 
(59) 
(60) 
(61) 



Here we have defined the delta- function 5 which incorporates a factor 1/ab. Based on these equations we can now 
justify the vanishing of the 55 and 05 components in the brane stress-energy (^5|). Such components, if non- vanishing, 
would appear on the right-hand sides of (^9|) and (|6(]) multiplied with delta- functions. We should, therefore, have 
corresponding delta-function terms on the left-hand sides of these equations. Since the scale factors a and b are 
assumed to be continuous, delta functions can only appear from second derivatives of these quantities with respect 
to y. However, there are no such terms in eq. (^9|) and (^0|). Hence, we conclude that the 55 and 05 components in 
eq. (^5|) must vanish. 

For the perturbations, we find to linear order 
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(67) 



Given those results, we can now return to the question of why the perturbations of the brane stress-energy tensors 
must have the specific form (57). We recall that all quantities in our metric ( f47| ) are even except the off-diagonal 
perturbation W which is odd under the Zi symmetry. From our continuity assumptions, delta- function terms in the 
perturbed Einstein tensor can, therefore, arise from first derivatives of W with respect to y and second derivatives 
with respect to y of all other quantities. Inspection of the above equations shows that such terms are absent in the 5i 
and 55 components of the perturbed Einstein tensor. Consequently, the corresponding components in the perturbed 
brane stress-energy should vanish. Furthermore, all terms in SG'j that could potentially lead to delta-functions are 
proportional to 8 l j. This implies that the anisotropic stress on the brane, which would contribute to the traceless 
part of eq. (p7|), must vanish. As a result, the traceless part of eq. (|67|) 



r)' 



-a 2 o\\ 3 



(68) 



only involves the bulk anisotropic stress as a source term. If the bulk anisotropic stress vanishes as well, as, for 
example, is the case for a perfect fluid in the bulk, one concludes that 



V>-0 + r = o 



(69) 



The quantity ip — <f> + T is the analog of the four-dimensional quantity -04 — (f>4 that also vanishes in the absence of 
anisotropic stress. This correspondence will be made more explicit in section V. 



IV. DENSITY FLUCTUATIONS ON THE BRANE 



A systematic study of density fluctuation in five dimensions requires solving the full set of five-dimensional equations 
of motion given in the previous section. However, for specific questions it might be useful to have some information 
about the metric restricted to the brane. For example, it is this restricted metric that is felt by matter which is 
confined to the brane. In this section, we are going to derive such equations on the brane starting from the general 
equations of motion above. 

For a Z2 even field the meaning of its value on the brane is quite clear. A Z2 odd field may jump across the brane 
so its value may have a sign ambiguity. To simplify the notation, we define the value of an odd field on the brane as 
the one that is approached from within the interval y S [0, R]. This is precisely the boundary value of the field as 
viewed in the boundary picture and, at the same time represents one half of the jump at the fix point. We also recall 
that the scale factors a, b and the perturbations 4>, tp an d T are even and hence continuous while the perturbation W 
is odd and may jump across the fix points. 
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Let us start with the background equations of motion (l^JT^]. We have already mentioned that the delta- function 
sources in eq. (|58|) and (|6l| ) have to be matched by the terms containing second y derivatives of the scale factors a 
and b. This leads to 



a o 



b' 



P 



(n) 



(70) 



These conditions, as well as the following ones, are valid at the brane positions y — y n where the upper (lower) sign 
holds for the brane n — 1 (n — 2). While the two other non- vanishing equations of motion do not contain delta 
functions, they can still be restricted to the brane. From the 05 component (|60|) we find 



P 



(n) = _ 3 « 



pin) + p (n) 



=F 2abr . 



(71) 



which represents to an energy conservation equation on the brane. Note, however, that, in addition to intrinsic brane 
quantities, this equation also involves the off-diagonal bulk stress-energy component r. This reflects the simple fact 
that the branes are not isolated systems but can exchange energy with the bulk. Finally, we should consider the 55 
component ( |59| ) . Restricted to the branes it results in an equation of motion for the values of the scale factors a and 
b on the brane given by 



a ab , , , a b 

+ kb 2 = 

a ab 3 



•3/ 



(72) 



An analogous procedure can now be applied to the perturbed equations. Observe that only the components SG° , 
SG l i , SG° 5 and 6G° i contain explicit delta- function terms. They should be matched by terms containing first y 
derivatives of W and second y derivatives of all other quantities. This leads to 

V/ = —W ± -ab (dp™ - Tp (n A (73) 
a 6 V / 

L Pi \ 1 1 . 

- Tp^j (74) 

W = [p (n> + P (n> ) v (n) (75) 

§ r (n) 

w - 7h ( 76 ) 

Interestingly, the last equation implies that the component (5T(™) 5 o of the brane stress-energy perturbation vanishes, 
as can be seen by comparison with eq. d57j). The component ST^ ^ = —5r^ n \ however, is non-zero and is, from 
eq. ([75[), (76) determined by 

SrM =T^P {n) (p {n) + P {n) )v^ . (77) 




We have, therefore, found an important additional constraint on the perturbed brane energy-momentum tensor (|57|). 
The quantity Sr^ is, in fact, uniquely fixed by the other components. For vacuum energy with = — p*-™' on the 
branes Sr^ n ' is zero, but it is generally non- vanishing otherwise. This is, perhaps, somewhat surprising since one could 
have expected that a purely four-dimensional stress-energy tensor on the brane (with all 5 components vanishing) 
should be allowed. Here we see that this is generally not the case. 

Next, we deal with the odd components 05 and 5i of the perturbed equations of motion given in ( |64| ) and j6^). 
Restriction to the branes leads us, after some algebra, to 

= - (l + wW) (t^l'n - 3^) - 3^ w^ n A tfW 

-2(1 + w {n) )v {n) a 2 (p + q) T tab (v + 20 - 6 {n) + — ^ -Lr (78) 

\ r J p^ n > 

and 
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v (n) |i 
6 2 
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(79) 



where w^ n > — p( ") /p(") and <5<") = 6pW denotes the energy contrast on the branes. These equations represent 
the conservation of energy and momentum for scalar perturbations, including possible energy-momentum flow from 
the bulk onto the brane or from the brane into the bulk. They should be compared with the corresponding equations 
in four dimensions, eq. (30) in [|TJ, taking into account that the variable 6 of (5lj is related to the peculiar velocity v 
via 9 = Eq. ( |78| ) differs from the four-dimensional result by the last two terms. They describe the energy flux 

from the bulk onto the brane. In eq. ( |79|) we have two new terms with respect to the four-dimensional equation. They 
describe momentum flux between bulk and brane. The coupling between bulk gravity and brane matter expressed 
via the above equations is one of the main results of this paper. It shows that, when considering scalar metric 
perturbations on the branes, the branes cannot simply be viewed as an isolated system but have to be considered 
together with the bulk environment. Practically, this implies that frequently one cannot simply copy four-dimensional 
formulae when dealing with physics on a brane that is embedded in a higher-dimensional space. Finally, we restrict 
the 55 component of the equations of motion, eq. (|(33|), to the brane. We find the following evolution equation for the 
perturbations projected onto the branes: 



& 2 (2V>-0)',.-3V>-3-0 + 3 --3- 
4 a \ o a 



ijj + 6kb 2 (ij + (t))+a 2 b 2 p {n ^ 



1(1 + 3wW) 



6q + 2ql If 3 (n) \ w ± ar (n) {n) 



(80) 



V. MATCHING TO THE FOUR-DIMENSIONAL EFFECTIVE THEORY 



In the previous subsection, we have derived a set of equations for the metric on the branes, essentially by restricting 
the five-dimensional equations of motion. These results may, for example, be useful to analyze the evolution of matter 
that is confined to the brane. However, the most important task is to extract predictions for structure in the late 
universe from our formalism of metric perturbations in brane-world theories. In this section, we will explain the basic 
steps in this direction. 

First, we should introduce the four-dimensional effective theory, describing physics at low energy, that is associated 
to our five-dimensional brane-world theory ( |53"| ) . It is this four-dimensional theory that described the evolution of the 
universe "today" and that is used for the interpretation of observational results. Theoretical predictions, originating 
from our brane-world theory, should therefore be formulated in terms of this effective theory. The five- and the four- 
dimensional effective theory are related by a vacuum state that constitutes a specific solution of the five-dimensional 
theory and should respect the symmetries that we expect to find in the four-dimensional theory. Specifically, four- 
dimensional Lorentz invariance implies that the vacuum metric should have the structure 

ds 2 = A 2 {y)dx ti dx v i 1 ^ ~ B 2 (y)dy 2 . (81) 

The functions A and B have to be such that this metric solves the five-dimensional theory in the vacuum configuration. 
For our five-dimensional theory (p3|), the simplest possibility is to have no stress energy in the vacuum which results 
in a flat vacuum metric, A, B — const. In five-dimensional heterotic M-theory the vacuum configuration is determined 
by certain potentials in the bulk and on the branes that involve the dilaton ifH^- In this case, A and B are non- 
trivial functions of y and the deviation from the flat vacuum metric is determined by the size of the so called strong 
coupling expansion parameter. The vacua proposed in ref. |^l],^2j are based on a vacuum configuration with specific 
cosmological constants in the bulk and on the branes and result in an exponential function for A in the coordinate 
system where B = const. Each one of these different vacuum states is associated with its specific low-energy theory. 
For the sake of simplicity and concreteness, we will here focus on the first possibility, namely the flat vacuum. This 
choice represents, at the same time, a good approximation for five-dimensional heterotic M-theory in the case of 
a small strong-coupling expansion parameter. The four-dimensional effective theory describes the dynamics of the 
collective excitations of the vacuum state. In our case, these excitation 3,16 given by £l four-dimensional metric c/4fj l v 
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and the modulus (3 describing the size of the fifth dimension. The vacuum metric with these collective modes put in 
has the structure 



ds 2 = e-Pg^dx^dx" - e 2/3 dy 2 , (82) 

where g^ and (3 are functions of x^. As usual, the effective four-dimensional description is valid as long as these 
functions are varying slowly enough. This is the case if all four-dimensional momenta are much smaller than the 
mass of the first Kaluza-Klein excitation around the vacuum state. In our case this mass is given by e~@ /2R. Let 



us, therefore, consider a five-dimensional evolution that is approaching the vacuum state (82). Even though the 
five-dimensional metric is then close to the vacuum metric it will still have small Kaluza-Klein excitations that can 
be described in linear perturbation theory. A useful way to extract the zero modes from such a five-dimensional 
metric with small contributions from Kaluza-Klein modes is to perform an average over the fifth dimension. Doing 
this systematically leads to the following four-dimensional effective theory associated to the brane- world theory ( |53| ) 
and the vacuum state (j32"l): 

R^u - \g^Ri = | (d^d v f3 - ^g^d(3 2 ^j + T ifiv (83) 
Vlf3 = ,h . (84) 



Four- and five-dimensional quantities are related by 



e 2 ? 


= - < 355 > 






= e 13 < g^ > 






-<T 11V > + 2Re0 


2 




71=1 


Ji 


= \e-^ <T 55 >- 





(85) 
(86) 

(87) 



Here, < • > denotes the average over the fifth dimension. As stated, this four-dimensional theory is a good description 
as long as all momenta are small compared to e~@ /2R, the mass of the first Kaluza-Klein mode. The Kaluza-Klein 
modes have decoupled from the above equations at linear order. However, due to the presence of the branes, the 
Kaluza-Klein modes cannot strictly be set to zero but have to be integrated out |^2|,^3|. This leads to higher order 
corrections to the above four-dimensional equations that are suppressed by powers of the four-dimensional Planck 
scale and that we have neglected. To be consistent with this approximation, the average < • > that projects out the 
Kaluza-Klein excitations should be considered meaningful only at the linearized level in Kaluza-Klein excitations. 

We would now like to apply the above general correspondence to our formalism for metric fluctuations. To do this, 
we need to assume a five-dimensional background solution that, for late time, approaches the vacuum configuration. 
Formulated in a four-dimensional language, this is the case if the four-dimensional Hubble parameter Hi = 04/04 
and $ are small compared to /2R. Furthermore, it is helpful to assume that the average < 700 > approaches 
one in this limit. This can always be achieved by a redefinition of time. We would like to explicitly work out the 
correspondence for the scalar sector in longitudinal gauge that we have focused on in this paper. The generalization 
to include vector and tensor perturbations is st raig htforward. Concretely, we apply the correspondence (pq)-(pq) to 



the five-dimensional quantities ([47[), (J54j), (|55j), (56) and ( |57| ) matching onto the four-dimensional quantities specified 
in (18), Jl9| ) and (pp[). Furthermore, we need to decompose the four-dimensional modulus /3 as 



P = x - r 4 ■ (89) 

where \ = x(t) is the time-dependent background and T4 = r^t, x l ) is the perturbation. The matching of background 
quantities leads to 

e 2x = < a 2 b 2 > 
a\ = e x < a 2 > 

pi = e~ x < p > + 



Pi = e x < p > 



2Re 2x 



2Re 2x 





(90) 




(91) 


n=l 


(92) 






71=1 


(93) 
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It is interesting to explicitly compute the background current J4 in the equation of motion ( |84| ) for the modulus x- 
It is given by 

Ji = 3 [pa - 3p 4 + 2e- x < <? >] . (94) 

With the above expression for p± and p 4 , this implies that the modulus \ has a runaway potential leading to a 
growing size of the fifth dimension. Therefore, the theory, as stands, will not stay in the range of validity of the 
four-dimensional effective theory. As is well-known, it needs further stabilization of the modulus (3 by means of a 
potential. In the context of string- or M-theory one expects such a potential to be generated by non-perturbative 
effects. 

The correspondence for the scalar perturbations reads 



r 4 = < r > (95) 
04 = < - T/2 > (96) 

^4 = < V> + r/2 > (97) 

1 2 

6 P4 = e-* <Sp + T P > +^-^ ]T(<5p (n) + 2lV")) (98) 

n— 1 

1 2 

6p 4 = e~* < Sp + Tp > +^-^ 5>P (n) + 21> (n) ) (99) 

n— 1 

cr 4 = e~ x < o > (100) 

„-3x„4 [ 1 2 

v, = -— ± <{p+p)v>+— T3 Y,(p {n) +P [n) > in) • (101) 
P4 + P4 2i?e' J ^ 

L n— 1 

In particular, we conclude that 

- 04 =< v- + r > (102) 



Hence, the four- and fivc-dimcnsional quantities that measure the presence of anisotropic stress are in direct corre- 
spondence to one another as they should. 

VI. CONCLUSION 

In this paper we have laid down a gauge-invariant formalism to describe metric fluctuations in brane- world theories. 
This formalism is a straightforward generalization of the well known formalism in four dimensions. It categorizes the 
perturbations according to their tensor properties with respect to the usual three-dimensional maximally symmetric 
space rather than a higher-dimensional space as one might have expected. This is a direct consequence of the brane- 
world nature of the theory which generally leads to cosmological backgrounds that are inhomogeneous in the additional 
dimensions. We have introduced a generalized longitudinal gauge in order to further study scalar perturbations. In 
the case of a five-dimensional model on the orbifold S 1 /Z 2 , on which we have focused, we have identified four scalar 
metric perturbations <p, ip, T and W. This has to be contrasted to the four-dimensional case where one only has 
two such perturbations. We have presented the evolution equations for these scalar perturbations which, mainly 
due to the dependence of the background on the additional dimension, are significantly more complicated than the 
corresponding four-dimensional equations. It is those additional terms, related to the non-linearity of the background 
in the additional coordinates, that encode possible new and interesting information about the formation and evolution 
of perturbations. Furthermore, given the gauge choices and assumptions about the continuity of the metric that we 
have made, we have determined the resulting most general form of the stress-energy on the brane. In particular, we 
have found that the perturbed brane stress-energy has to have vanishing anisotropic stress and that its 05 component 
is non-zero. We have applied our formalism to calculate the matching conditions (Israel conditions) for the five- 
dimensional metric restricted to the branes. Among other results we have derived the equations describing energy 
and momentum conservation for metric perturbations on the brane. As is expected on physical grounds, they illustrate 
that the brane cannot be viewed as an isolated object but is subject to energy and momentum flow between the bulk 
and the brane. Finally, we have shown how the five-dimensional formalism for metric fluctuations can be matched to 
the known four-dimensional one in the limit where the brane- world theory has an effective four-dimensional description. 
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This allows one to extract predictions for structure in the late universe originating from brane-world theories. We 
hope to address this problem more explicitly in a future publication. 

As this manuscript was prepared for submission, ref. fl54| , ]55| and |5(| appeared which have some overlap with the 
present paper. Cosmological perturbations are also discussed in 
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